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ABSTRACT 


This report is concerned with theoretical predictions 

of plume interference {signal attenuation) for the Space 

1 

Shuttle Range Safety system by the Solid Rocket Booster 
(SRB) exhaust plvunes. The signal propagation is calculated 
using a split operator technique based upon the Fresnel- 
Kirchoff integral, using fast Fourier transforms (FFT's) 
to evaluate the convolution and treating the plume as a 
series of absorbing and phase- changing screens. Talanov's 
lens transformation is applied to reduce aliasing problems 
caused by ray divergence. 

Single-plume flow fields supplied by Lockheed Missile 
and Space Company were used as inputs. Three-dimensional 
effects at low altitudes were simulated by superposing the 
electron density and collision frequency of two independent - 
plumes; at higher altitude the plumes are treated as 
merged plumes from an equivalent nozzle. Simulated plume 
turublence was included by constructing realizations of the 
displacement field in which possible correlations with 
axial distance were taken into account. 

Comparisons of calculations of the mean flow plume- 
induced loss with field data from Pro;]ect See-Thru show that 
the predicted signal loss for Titan III-C agrees fairly well with 
observations. The effects of simulated turbulence upon the 
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mean attenuatxon are small . Calculations for the Space 
Shuttle SRB plume predict attenuation levels exceeding 
60 db at tail aspect at altitudes above about 20 km. The 
attenuation pattern becomes broader and deeper at higher 
altitudes, with maximum signal loss at the highest altitude. 
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1.0 INTRODUCTION 


In its capacity as the operating agency for the 
Eastern Test Range (ETR) , the United States Air Force 
requires that launch vehicles and missiles fired from Cape 
Canaveral be equipped with a Range Safety System. This 
system should be capable of destroying the vehicle in flight 
on command from the ground when, in the judgment of the 
Range Safety Officer, life or property ad^jacent to the 
test range is threatened by malfunction or course deviations 
by the test vehicle. The Space Shuttle, whose first launch 
is scheduled for 1979, will carry a Range Safety system, 
and the George C. Marshall Space Flight Center, NASA, has 
been charged with the design and qualification of the system. 

The SRM's burn a high-performance aluminized composite 
propellant. Previous experience with other large vehicles 
of this type (Minuteman, Polaris/Poseidon, Titan III) 
indicates that substantial signal loss and plume-induced 
noise are experienced when the line-of -sight signal passes 
through the exhaust plume . The Shuttle SRM's are the largest 
solid-fuel motors ever to be carried by a U.S. launch 
vehicle, and the exhaust gas is expected to be highly ionized 
in relation to most liquid-propellant plumes. A potentially 
serious signal transmission problem could be present when 
the plume intercepts the signal path. 
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The study reported here was undertaken to predict 

plume effects on the Range Safety signal. The last time 

that a similar problem was addressed for large rockets was 

about ten years ago when the Saturn and Titan III were under 
* 

development. In the meantime, computational procedures 
both for exhaust plume flow fields and chemistry and for 
forward electromagnetic wave propagation in lossy media have 
undergone considerable advances. These new procedures have 
been applied to our predictions. This study can therefore 
be considered both as an engineering study of an operational 
problem and as an advance in the state of the art of plvime 
signal interference analysis. 


*Development and analysis has continued for tactical missile 
systems, much of it at the Naval Weapons Center. 
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2.0 THE SPACE SHUTTLE RANGE SAFETY SYSTEM* 

V 

The Space Shuttle Range Safety System is simxlar xn 

concept and operation to that previously employed on the 

Saturn V in the Apollo program. The command signal is a 

frequency modulated, frequency-shift keyed (FM-FSK) sxgnal 

with a nominal central frequency of 416.5 MHz. The signal 

consists of a series of 11 tone pairs. The first 9 tone 

pairs alert the system, and the last 2 carry the signal 

Itself. The tone pattern is shown schematically in Figure 1. 

The frequencies used range from 7.35 KHz to 13.65 KHz at 

intervals of 1.05 KHz. These frequencies are well outside 

the, bandwidth of previously determined plume noise, so that 

plume noise should not have a significant effect on the 
** 

signal. 

MoTonting of the Range Safety antenna is shovm in 
Figure 2. Each SRM carries two antennas mounted 180° apart 
and as far forward as possible (about 40 meters forward of 
the engine exit plane) . Two additional antennas are mounted 


*Most of the information in this section was furnished by 
Lee Malone, NASA/MSFC. 

* 'k ' 

Plume noise can have an effect on the S-band communications 
link over which some of the information required to make 
a range safety decision is transmitted. We have not 
addressed this problem in any depth in this study, although 
similar procedures should be applicable. H.A. Poehler has 
pointed out that flame noise can broaden the Range Safety 
Signal spectrum, possibly to the point where the spectrum 
could overlap with the result of erroneous detection. 
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Figure 1. Space Shuttle Range Safety Signal Modulation Envelope 












on the external tank, also well forward. Coupling with 
vehicle skin currents is not known, but is probably fairly 

* ’ j / ‘ 

strong for this complicated configuration. The hear . field 
antenna pattern has not been measured, although far field 
patterns are known. In this work we have assumed 
the antenna pattern to be isotropic. 
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3,0 PLUME-SIGNAL INTERFERENCE PHENOMENOLOGY 


In this section we address the signal interference 
problem from two different points of view. We first adopt 
an empirical approach in which we review selected portions 
of the available data base, looking for trends (or lack 
of trends) with altitude and signal frequency. We then 
switch to an analytic approach, considering factors 
affecting RF propagation in a non-uniform, lossy dielectric 
medium, the relation of plume properties to propellant 
properties and operating conditions, and determine those 
effects which are theoretically expected to determine the 
observed signal levels. These considerations determine the 
computational approach adopted in the following section. 

3 . 1 The Empirical Data Base 

Because of the importance of plume interference to 
communications and telemetry, considerable attention has been 
given to it and a substantial body of data exists. Data 
for large vehicles are available from both Eastern and 
Western test ranges (ETR and WTR) . Victor (1975) has 
summarized the results of many model tests and small vehicle 
flight tests. Poehler (1969 a, b, c) has summarized the 
best available ETR measurements, covering studies performed 
over the period 1963-1967. Because many of the older data 
contain unknown signal loss contributions from antenna 
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Figure 3. Composite representation of variations of P-band 
signal strength for Titan III-C launch in 
Project See-Thru. 


pattern "holes" and , pointing errors, they are generally less 
reliable than the later data. 

The most complete body of data for a single vehicle is 
that associated with Project See-Thru -(Poehler, 1966, 1967, 1968) , 
a study of plume-signal interference effects associated 
with Titan III-C launches. These data cover several aspect 
angles with respect to the vehicle roll axis as well as 
several signal frequencies. In addition, the morphology 
of the Titan III-C 0^^ stage plume is somewhat similar 
to that of the Space Shuttle SRM plumes, since in each 
case the two solid motors are parallel-staged with the 
upper stage. The effects found in these measurements are 
representative of those shown by plumes from other large 
solid-fueled vehicles, and the quoted attenuation levels are 
probably more quantitative. 

Figure 3 is a composite representation (Poehler, 1966) 
of the sequence of events observed at P-band (250 MHz) , in 
which observations from several fixed ground locations have 
been combined. (Few of the actual data plots are of this 
quality, and none will show all the effects shown in this 
composite.) No signal effects attributable to plume inter- 
ference are observed during 0^^ stage burn until about 75-90 
seconds after liftoff. The aspect angle is high during 
this period, so that the line of sight does not intersect 
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the plume. At different times for each ground site, the 
line of sight intersects the plume and the signal becomes 
noisy and attenuated. Plume noise and plume attenuation 
are coupled; a noise-free signal loss is due to antenna 
pattern, pointing errors, or other phenomena. The aspect 
angle at a given ground site may either increase or decrease 
during the period from 80-110 seconds. The attenuation 
level typically increases with decreasing aspect, suggesting 
that the longer the signal path through the plume, the more 
greatly it is attenuated. First stage ignition (110 seconds) 
causes increased attenuation and SRII separation (122 seconds) 
produces a drastic drop in signal level. Additional periods 
of attenuation are observed during 1st and 2nd stage burn 
and at 2nd stage ignition. 

Data were obtained at S-band (2.22 GHz), C-band 
(5 GHz) and X-band (9.1 GHz), as well as at P-band (250 MHz). 
Because the polarization characteristics of the X-band 
antenna v;ere not determined, the X-band data are reported 
as "signal loss" rather than as "plume attenuation". When 
one attempts to interpret these data in terms of dependencies 
on various parameters (aspect, altitude, signal frequency) 
several complicating factors arise. Because measurements 
were taken at fixed ground locations, variation in signal 
characteristics at any given site contain the combined 
effects of aspect angle and altitude. Only by comparing 


10 



results from different sites xs it possible {in principle) 
to infer aspect angle dependence at a given altitude- In 
addition, measurement problems at times intrude; for 
example, some stations show unexplained decreases in signal 
levels below the predicted free-space curve, and rapid 
fluctuations of up to 20 db in signal level occur during 
periods of plume interference. For these reasons Poehler 
has reported many of his data in the form of upper or lower 
bounds to the actual signal level. 

Summaries of the attenuation measurements from Project 
See-Thru are given in Tables 3-1 and 3-2. The results in 
Table 3-1 are Poehler 's summary of data from ETR Test 
8275/2250; Table 3-2 is our interpretation of the P-band 
data from Tests 3656/6025 and 6020/6546, and could be in 
error by more than 5 db. 

These data have been examined for trends in altitude 
and signal frequency. In Figures 4 and 5 the data 
obtained near 80 seconds and 110 seconds after liftoff 
(altitudes near 20 and 36 km, respectively) are plotted as 
attenuation versus aspect angle. All frequencies are shown. 
While at each altitude a trend with aspect angle is evident, 
there is no clear separation of data in terms of signal 
frequency (except at X-band, where the calibration contains 
an unknown correction for antenna polarization losses) . 
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Table 3-1 




X-Band Signal Loss**, Zero Stage 


Aspect Angies* 

Signal Loss 

Test Time 

Site 

a = 180 -0 

0 




> 15 

70-150 

0 db 

0-110 sec 

UC-2 

16 

339 

> 10 

+ 65 

19.1 

10-8 

245-223 

< 30 

95-110 

UC-11 

9 

186 

30 

108-110 

19.1 

4-5 

199-188 

40 

100-110 

UC-13 



TABLE VI 





C “Band Attenuation, Zero Stage 




Attenuation 

Test Time 

Site 

> 15° 

50-150 

0 db 

0-110 sec 

GIF 

11 

346 

> 10 

+ 74 

UC-13 

6 

188 

> 20 

+ 110 

UC-13 



TABLE Vn 





S-Band Attenuation, Zero Stage 


> 25° 

100-137 

0 db 

0-110 sec 

STS 

> 16 

70-150 

0 

0-110 

UC-2 

> 15 

50-150 

0 

0-110 

GIF 

15 

306 

> 8 

+ 75 

UC-11 

10 

346 

> 12 

+ 74 

UC-13 

8 

223 

> 10 

+ 110 

UC-11 

5 

188 

> 20 

+ 110 

UC-13 



TABLE Vm 




P 

-Band Flame Attenuation, 

Zero Stage 


V 

o 

329-306 

0 db 

0- 75 sec 

UC-11 (and UC-13) 

8 

223 

15 

+ 110 

UC-11 

5 

191 

15 

+ 105 

UC-13 

2 

231 

20 

+ 90 

UC-13 
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Table 3-2 


P-Band Telemetry 

Data, 

Titan III-C, 

AFETR 

Vehicle Time 

Atten 

Aspect 

Site 

3656/6025 0-110 

none 

>18° 

Tel 2 

105 

10 

15° 

UC-12 

80 

18 

4° 

UC-18 

95 

~10 

9° 

UC-18 

110 

15 

11.5° 

UC-18 

6020/6546 110 

~5 

14° 

Tel 2 

80 

-5 

13° 

UC-7 

95 

~10 

13° 

UC-7 

110 

>18 

11.5° 

UC-7 

80 

<16 

7° 

UC-18 

95 

<18 

8° 

UC-18 

110 

<22 

8° 

UC-18 
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Figure 4. Attenuation levels as a function of polar aspect 
angle for Titan III-C near 80 sec. TALO 
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Figure 5. Attenuation levels as a function of aspect 
angle for Titan III-C near 110 sec. TALO. 



Figure 6. Attenuation levels at P-band {-^250 MHz) as a 
function of polar aspect angle for Titan III-C 
at various altitudes. 


TALO 


Figure 6 shows the P-band data grouped into three altitude 
regimes; the attentuation here appears to be marginally 
greater at higher altitudes at a given aspect, but again 
no clear trend is evident. 

At this point we shall anticipate the results of 
the following section and state that in general we expect, 
on theoretical grounds, an effect of both altitude and signal 
frequency upon the observed signal level. That there is no 
clear-cut dependence in the best available data is on its 
face somewhat perplexing. We conclude that a careful 
accounting of the various factors at work -is required, since 
opposing trends and subtle interplays among them may be at 
work . 


3 . 2 Theoretical Considerations 

Signal transmission through an exhaust plume is an 
example of electromagnetic wave propagation through a 
lossy dielectric medium. It is convenient to divide 

the discussion into separate sections on the properties of 
the medium (the plume) and its effects upon a propagating 
wave. 

The plume interferes with an RF signal because it 
contains free electrons and ions. These are produced by 
thermal ionization of alkali metal impurities in the pro- 

j 

pellant. These levels tend to be high (50-100 ppm) in 
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composite propellants because of the chemical similarity be- 
tween sodium and potassium and the ammonium group present 
in the oxidizer (ammonium perchlorate) , which makes puri- 
fication difficult. Ionization occurs both in the combus- 
tion chamber and in the hot portions of the external flow. 
Once produced, electrons are removed by recombination with 
ions and by attachment to neutral atoms (principally 
oxygen and chlorine) . The dielectric constant e depends 
upon the propagation frequency determined by the electron 
number density and their collision frequency with the 

neutral species (the fractional ionization level is always 
small, as is usually the ion contribution to e ) . An 
explicit relation is 


2 

w 



U^v/o) 

-I- -1 tl 

2 ^ 2 , 2 
V W + V 


(3-1) 


where is the plasma frequency 

47t e N q ^2 

(0 =; -^=3.2x10^ N (sec if N in cm ) 

P e e 


(3-2) 


( e is the electron’s charge, and its mass) , and 
0 ) IS the circular frequency of the transmitted radia- 
tion. 

A more directly useful quantity is the refractive 
index, n = /z . Suppose that a wave travels a small 
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distance 6s xn a medxum whose refractive index is n. The 
ratio of the electrxc field at the end of this path to its 
value for free-space propagation over the same distance xs 

E^/E^ = exp{x k(n-l) 6s} (3-3) 

where k = w/c is the wavenumber. Because e and thus n 
is complex, the exponent contains a real and an xmagxnary 
part: 

E^/E^ = exp {i k (n^-l) 6s - k n^^ 6s) (3-4) 

The imaginary part of the exponent represents a phase shxft 
of the fxeld which results from the change in wave phase 
velocxty as the refractive index varies from unity. The 
real part is negative and represents an attenuatxon of the 
wave, or in other words, a transfer of energy from the field 
to the medium. Both the real and imaginary parts of the 
refractive xndex increase wxth increasxng electron densxty 
and decrease with xncreasing signal frequency. The depen- 
dence on collisxon frequency is more complxcated. Explicit 
expressions for the real and imaginary parts of n are 

1 - X t + r (1-X) ^ 

2 (1 + Z^) I 4 (1 + Z^) 
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and 



where 

X = (3-7) 

and 

Z = v/w . (3-8) 

The parametric variation of n and n. are shown in 

r 1 

Figure 7 (Kelso, 1964) . 

The electron density and collision frequency, and 
thus the refractive index, vary with position and time 
within the exhaust plume. This variation is determined by 
the gas dynamics of the exhaust gases as they expand to 
ambient pressure, mixing of the exhaust with the atmosphere, 
and chemical reactions. The basic structure near the 
vehicle is determined by inertial and pressure forces as 
the exhaust expands on leaving the rocket nozzle, and for a 
single plume takes the form of a core expansion surrounded 
by a shock layer. Superimposed on this structure is a 
region where turbulent mixing of the exhaust and the ambient 
atmosphere takes place. This mixing region grows as the 
gases move away from the vehicle, while the pressure forces 
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Figure 7 - Real part ]x and imaginary part X of the 

complex index of refraction as a function 
of the parameter X = 03p2/u2 for several 
values of the parameter Z = v/o). (Figures 
from Kelso, Radio Ray Propagation in the 
Ionosphere . ) 
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dimxnish. Far downstream the exhaust gases are at the same 
pressure as the ambient atmosphere, and the further develop- 
ment of the flow is controlled by turbulent mixing. 

f 

Within the mixing region chemical reactions can occur 
which release heat to the flow field and cause additional 
ionization. Because the exhaust gases are generally fuel- 
rich, for reasons associated with motor performance, 
secondary combustion or afterburning can occur as they mix 
with air. The chemistry is fairly complex, but reasonably 
w'ell understood. The reactions are generally out of equilibrium, 
so that ionization levels are determined by a complex set 
of rate processes. 

Computer codes have been developed which treat the 
combined effects of expansion gas dynamics, mixing, and 
finite-rate chemistry to evaluate time-averaged electron 
density and collision frequency as a function of position in 
the plume. Some results of calculations performed at Lock- 
heed Missiles and Space Company in support of this study are 
shown in Figures 9 and 10. At higher altitude, the plume 
is broader, longer, and less dense. For .typical antenna 
locations the aspect angle at which the plume first inter- 
cepts a ray from the antenna increases with altitude. At 
the same time the difference of the refractive index from 
unity tends to decrease. Qualitatively, we exnect the plume 
to be large but more transparent at higher altitude. Which 
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effect predominates -will depend upon vehicle size, propel- 
lant, trajectory, and signal frequency. 

Detailed calculations (e.g. , Thomson, 1962) indicate 
that in the case of large and relatively highly ionized plumes, 
such as those from the SRM's for Titan III and the Space Shuttle, 
rays which enter the plume are both strongly absorbed and 
appreciably refracted. As a result, the plume casts 
a very dark shadow on the transmitted signal. Diffraction 
of rays which pass a short distance outside the shadowed 
region, though generally a weak effect because of the 
large value of the ratio of a characteristic plume radius to 
the signal wavelength, can produce signal levels in the 
shadowed regions which are several orders of magnitude greater 
than those associated with direct rays (Figure 10) . 

The interpretation of the signal near tail aspect, 
deep within the plume shadow, as a diffracted signal, has 
several implications as to the processes in the plume which 
determine the signal level. At any given signal frequency 
lov/ enough that the plume is essentially opaque, the signal 
in the shadowed region is a sensitive function of the 
effective plume diameter and of conditions near the edge of 
the plume. The edge region is one of high mean electron 
density gradient, a consequence of diffusion and of the 
chemical processes which remove electrons from the flow. 
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Figure 10 . C-band attenuation calculations for ELP 

propellant in flight at 25,000 ft. altitude. 
(Flight data are shown for the altitude 
range 15,000 to 31,000 ft.) Data from 
Victor (1975) . 
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It is also one in which the turbulent flow is intermittent. 

.1 ’■ 

Large', edjlaes'.^spo.radically convect highly ionized gas from 
the interior of the mixing zone into and out of this 
region; the actual plume edge will be at any instant an 
irregular structure dominated bv this large eddy structure. 
Figure 11, a shadowgraph (Laufer, 1975) of conditions in a lab- 
oratory mixing layer between fluids of different density, illus- 
trates the type of structure which we expect in the plume 
edge region. At a location where the mean electron density 
is (say) 10 , with interior densities of 10 , the mean 

density is more likely a consequence of the local condition 

10 * . .8 

being 10 one percent of the time than of its being 10 

all the time. The RP signal responds to the instantaneous 

structure of the plume. 

The interpretation of the shadow zone signal as a 
diffracted signal also leads to an opposite dependence on 
frequency and altitude to what one would have expected for 
a transmitted signal. At higher frequencies (yet low enough 
that the plume is opaque) the signal propagation is more 
nearly geometric, so that one expects a reduced tail-aspect 
signal as the frequency is increased. At higher altitudes 
(yet not so high that the plume becomes transparent) the 
shadowed zone is broader, and the filling in of the shadowed 
region by diffraction is reduced. On the other hand, to the 
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Fiqure 11. Shadowgraph of a mixing layer produced 
by helium (upper side) at = 1015 cm/ 
sec and nitrogen at U 2 = 384 cm/sec at 
a pressure of 8 atm (from Laufer) . 
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extent that ed^^e gradients are somewhat smooth, the plvime 
shadow. ALialtlbe somewhat less broad at higher frequencies. 
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4.0 COMPUTATIONAL PROCEDURES 


In order to calculate plume-induced signal attenuation 
we should employ a physical optics procedure which accounts 
for the combined effects of absorption, refraction, and 
diffraction. In this section we describe such a procedure. 
Similar techniques have previously been applied to RF 
propagation in the ionosphere (Chestnut, 1972; Watson and 
Thomson, 1972), laser beam propagation and thermal blooming 
in the lower atmosphere (Fleck, Morris, and Feit, 1975; 

Berger, et al. , 1977) and sound propagation in the ocean 
(Flatte and Tapper t, 1972) . The application to plume- 
signal interference is novel, and we therefore describe the 
procedure in some detail . The procedure is developed in a 
heuristic fashion; a more formal development .is given in 
the Appendix. We begin by discussing the wave propagation 
problem in free space. A treatment of absorption and re- 
fraction in a material medium is then introduced. Finally, 
we discuss procedures employed to reduce numerical problems 
arising from ray divergence and other factors . 

Before beginning, we may refer back to Figure 10 to 
note that the deflection angles of rays refracted by or 
diffracted around the plume are small (generally less than 
10-20°). Under these conditions propagation and scatter 
take place in a generally forward direction. We may therefore 
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characterize the propagating wave by a wave vector k which 
is in the nominal direction of propagation, and employ 
small-angle approximations (cose 1, sine a* tan 6 6) 

throughout. This behavior considerably simplifies the 
analysis and is appropriate so long as we consider a trans- 
mitted (or forward-scattered) signal rather than reflection 
or backscatter. 


4.1 A Free-Space Propagation Algorithm 

Kirchoff's integral relation, the mathematical 
statement of Huggens' principle, is the basis of physical 
optics diffraction theory. As originally developed it is a 
scalar theory, although vector formulations have been 'developed 
(Jackson, 1962) - Since plume dimensions are typically large 
with respect to a wavelength and vector wave corrections are 
important only v^hen v/avelength and aperture dimensions are 
comparable, the scalar formulation is adequate for our 
purposes. Let us assume that we know the field amplitude 
Ip on a plane Z = normal to the nominal propagation 
path, which we take to be along the Z axis. Then the field 
at points on a plane- Z = Z^ some distance further along 
the Z axis is given by the relation 


oo oo 




ikr -1 - n 
e r dx dy 
o -^o 


—CO — CO 


(4-1) 
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where r is the distance between ' 

Since r is 

r = ’ (4-2) 

the term multiplying the known field in (4-1) is, at a 
given distance between planes, a function of (x^^-x^) and 
(Yi-Yo) . Thus, the field at Z = can be expressed as 

a convolution, 

CO CO 

^(*1-^0' 5’i-yo> '^"‘o ‘*^0 • 

— CO —CO 

The easiest way to evaluate this convolution, when the number 
of points in each plane is large, is to multiply the Fourier 
transforms of and F: 


i)^{k^,ky) 


z=z 


1 




z=Zo 


♦F(k , 

X 



(4-4) 


where the tilde denotes a transformed function. Past Fourier 
transforms (FFT's) allow these operations to be done rapidly, 
with a number of steps proportional to log 2 (N^ Ny) 

instead of the operations required by direct 

integration at every point. 

An efficient computational procedure for evaluating 
forward wave propagation in free space thus consists of the 
following steps: 
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1 . 


Evaluate the Fourier transforms of 
and the propagation kernel F. 

2. Multiply the transforms at each point in 
k-space to obtain the transform of their 
convolution 

3. Evaluate the inverse transform to obtain 

^ (^I'^l^ * 

The process is repeated for however many planes at which one 
wants to have the field explicitly displayed. 

For the Kirchoff integral, in the form of Equation (4-1) , 
the Fourier transform of the propagation kernel F cannot 
be expressed in any simple analytic form. However, by 
invoking the Fresnel approximation 


F = 


ikr 

09 ■— 

r 2ir 


+ ik (z 


we immediately obtain 


exp + (yj^-Y^)^l/2 


(4-5) 


F = 


ik(z^-z^) 

exp |i[z-2^) (^^+kp/2 k 


} 


(4-6) 


Here we assume that the deviation from forward propagation 
IS small enough that we may neglect differences in [r| as 
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but do 


they affect the modulus of the V7ave reaching , 
not restrict these differences to be small with respect to 
a wavelength, so that phase effects are important. Since we 
are considering essentially forward propagation, this version 
of the small-angle approximation is adequate. 

A representative application of this procedure 
to simple diffraction problems is shown in Figure 12 
and 13. In the second of these we show the calculated 
far-field (Fraunhofer) diffraction pattern of a rectanular 
aperture, and in the first a photograph from Jenkins 
and White (1957) . 

4 . 2 Effects of a Dielectric Medium 

Instead of propagation in free space, we want to 
treat propagation in a medium whose refractive index varies 
with position. Let us first consider a thin slab of inter- 
vening material and see how we might modify the treatment 
for this case. Referring to Figure 14, we consider a 
source field on the plane z = z^ , a slab of dielectric 
material centered about the plane z = z^ , and ask for the 
resulting field amplitudes on the plane z = Z 2 - Beginning 
with the field at z^ , we calculate the field incident upon 
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DIFFRACTION PATTERNS OF RECTANGLES 



Figure 12. 


Fraunhofer diffraction of 
a rectangular aperture (from 
Jenkins and White, Fundamentals 
of Optics) . 



2X10 

Figure 13. Calculated Fraunhofer 

diffraction pattern of a 
rectangular aperture showing 
contours of constant intensity. 
Calculation is symmetric about 
the abscissa. 


ORIGINAL PAGE IS 
OF POOR QUALITY 




the slab by applying the free-space propagation algorithm 
between and . The wave emerging . from the slab, 

according to the discussion in Section 3.2, is a wave 
which IS corrected in amplitude and phase: 

Tp(x,y,z^ + h/z) = i|j(x,y,z^ - h/z) exp ^i [n^ (x,y) -l] k 6z 

- [n^{x,y) k' 6z] + i k 6z ^ . (4-7) 

The last term in the exponent is a position-independent phase 
correction and can be dropped. The second term represents 
absorption by the slab, and the first term a position-dependent 
phase shift. The phase fronts of the emerging wave will be 
tilted with respect to those of the incident wave to the 
extent that the real part of n varies with position; this 
IS equivalent to a change in direction of the wave, and 
corresponds to refraction by the gradients in the slab 
perpendicular to the nominal propagation direction (the z 
axis) . By applying the free-space propagation of the emergent, 
wave at , we can evaluate the field at Z 2 - 

The extension of this procedure to a thick medium is 
straightforward. We divide the medium into a number of 
slabs p separated by distance Az^ . Between each adjacent 
pair of slabs we allow the wave to propagate as if in free 
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space using FFT'^s to evaluate the Fresnel-Kirchof f integral. 

At each z station z . we apply an amplitude and phase 

• J 

correction at each of a preselected set of grid points, 
appropriate to the local refractive index and the slab 
thickness Az^ . The corrected wave is then propagated to 
the next location this process of alternate 

propagation and correction is completed until the entire 
medium has been traversed. Once through the medium, a 
final propagation step gives the transmitted field at any 
desired distance. The process is shown schematically in 
Figure 15. 

It can be shown (e.g., Fleck, et al . , 1975) that this 
algorithm is a solution, correct to second order in Az , 
of the parabolic wave equation 


2ik II = vj iji + k^(n^-l) ip 


(4-8) 


The assimiptions employed in deriving Equation (4-8) are 
equivalent to the Fresnel approximation. In operator form 


the solution is formally expressed as 

,2 

ik 


1 Az V 




4k 


Az(n^-l) - 


i Az V 
4k 


^3 


(4-9) 
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Figure 15. Schematic process of evaluation of 
signal propagation in an absorbing 
and refracting medium. 
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where 


(4-10) 



rj+1 


(n-1) dz . (4-11) 


The algorithm is also applicable to nonlinear wave propaga- 
tion, in which the local refractive index is a function of 
the local field amplitude, and has been widely used in that 
field. 

At this point we have available a computational 
procedure which includes the effects of refraction, absorption, 
and diffraction in a single calculation- Because each 
point on the plane z = affected by every point on 

the previous plane z = Zj , the procedure includes multi- 
path effects. The phase-correction step is applied at 

^y Poi^^ts in each plane, while the convolution operation 
requires a number of steps proportional to log 2 (N^Ny) . 

For a large number of points, the total computation time 
is therefore roughly proportional to , and the 

field amplitude is determined at each of N N points in 

X y 

the target (or ground) plane. 
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4.3 Some Numerical Detaxls 


Several problems arise in applying this procedure to 

propagation in an exhaust plume to or from a wide-angle 
* 

antenna. These are primarily associated with the divergence 
of the rays from the antenna. The convolution of Equation 
(4-1) in principle requires an integration over the entire 
plane z = . In approximating it by a finite-range 

discrete Fourier transform, we leave ourselves open to errors 
arising from the finite mesh spacing in both physical (x) 
and Fourier (k) space. These mesh sizes are related as 


6k 


2 IT 

N 6x 

X 


(4-12) 


When we limit the range in physical space, we eliminate 
the low-wavenumber Fourier components lying between 6k and 
k = 0 . This imposes an implicit assumption that the function 
is periodic in x space. A similar assumption of periodicity in 
k space results from the fact that we choose a finite mesh 
size in x space. Because the function is implicitly 
periodic, a result which runs into the mesh boundaries in 
either space results in a calculation in one periodic cell 
★ 

A well-known reciprocity theorem allows us to treat the 
propagation in whichever direction is convenient, since 
the signal losses along any path in either direction are 
equal . 
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containing contributions which are in reality associated 
with an adjacent cell. In Fourier space this behavior is 
called aliasing ; high-wavenumber components appear as 
effective low-wavenumber components. In physical space it 
IS called imaging ; effects which should appear at large x 
appear at small x. It is as if ghost or image solutions 
developed in the computational cell adjacent to the one in 
which we are interested, and our solution contained contri- 
butions from these images. Thus, once the rays from the 
antenna reach the edge of the grid, it is as if we were 
seeing contributions from antennas not only at (say) 

(x=0, y=0) but also from antennas at (x = N 6x, y=0) , 

(x=0, y = Ny 6y) , (x = 6x, y=0) , (x=0, y = -N^ 6y) , 

and so on (see Figure 16) . Brigham (19 74) gives a very 
clear description of these aliasing and imaging phenomena 
(Gibbs phenomena) and shows how they are related to the 
properties of discrete Fourier transforms. 

Since we have no a priori way of eliminating these 
image contributions to the computed field once they occur, 
we seek a procedure which minimizes or eliminates their 
occurrence- A coordinate transformation due originally to 
Talanov (1970) has been found useful in treating focused 
laser beams, where similar problems arise. Talanov 's trans- 
formation is readily applied here. We write the field on 
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Figure 16. Imaging (leakage) effects caused by discrete resolution on 
( a limited grid. 



some plane z = Zj as the produce of some function 
ij^'(x,yrZj) and a spherically symmetric phase factor: 


i'k(x^+y^)/2 z. 
»|j(x,y,Zj) = 4^' (x,y,Zj) e 


(4-13) 


where rs the focal length of the lens which produced 

the wave front represented by the phase factor (positive 
for converging beams and negative for diverging beams) . Then 
the freely propagating wave at z = can be expressed 

as 


4;(x,y,z ._^, ) = 


j+1 1-Az/z^ 


4 ;' 


X 


, z + 


Z . , T - Z 

1+1 3 


1-Az/z^ ' 1-Az/z^ ' 1 1-Az/Zj 

j . (4-14) 


X exp|ik (x^+y^) / [2 z^{l-Az/z^)] 


Applying this transformation involves the following steps; 

1. Divide the known field 4) at each point of the 
plane z = z^ by the phase factor to obtain 
the residual field 4^ ' . 

2. Calculate the field ip' at the plane z = 
using the algorithm of Section 4.1. 

3. Rescale the field amplitude and the mesh according 
to the relations 
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ijj" = Ij;' 


1-hz/z^ 


x' = x{l-Az/z^) 


y' = y (1-Az/z^) 


' 

> (4-15) 


4. Correct the phase of at each point of 

the plane z = by multiplying by the 

phase factor 

iJc(x'^+y’^)/2(z^-Az) . 

e 

The phase and amplitude corrections introduced by the medium 

are then introduced as in Section 4.2 and the calculation 

is continued to the next z station. 

This transformation converts a spherical wave (in the 

small-angle approximation) into a plane wave. The mesh 

on which the field is evaluated thus grows or shrinks as 

the calculation proceeds, as the focal point is behind or in 

front of the station being considered. The propagation algorithm, 

given a uniform plane wave at z = (which has only a k = 0 

Fourier component) maintains it as a uniform plane wave 

1 

at all successive stations. The combination of this algorithm 
and the Talanov transformation also maintains the character 
of a (small-angle) spherical wave. It will also preserve a 
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Fraunhofer diffraction pattern, provided that the calcula- 
tion IS initiated at som-- staiion at which the wavelength 
IS very sma3 ■* with respect to the product of a pattern 
dimension tir.es a diffraction angle. 

In applying this j'rocedure to the plume-signal inter- 
ference problem, we generally assume the focal point to be 
at the antenna. For most conditions of interest, the 
antenna is located far enough forward on the vehicle that 
the angle subtended by the plume is small (of order 20®), 
so that we can consider both the plume and the space for 
some distance around it within the confines of the small- 
angle approximation. Three-dimensional flow fields can be 
treated as easily as axially symmetric ones, and local in- 
homogeneities are accommodated to the resolution of the mesh. 

Since the outer regions of the plume are reflecting, 
an additional ray divergence is introduced over and above 
that of the initital spherical wave. Depending on the re- 
flection angle, these rays can sometimes get to the mesh 
boundary and cause image effects. It is usually easy to 
recognize these if they are not too strong, but they are an 
artifact of the calculation. Likewise, rapid variations in 
plume conditions in a plane introduce high-wavenumber 
Fourier components and the possibility of aliasing in Fourier 
space. It is usually desirable to introduce a rolloff 
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factor which reduces the contribution of the high-k components 
(and thus may eliminate some fine structure from the result) . 
An apparently satisfactory formulation is 

*corr = Wc + 

where k = tt/6x is the maximum wavenumber carried, A 
max 

similar reduction can be applied to the values of the field 
near the edge of the grid in x-space ; here the factor is 
applied to the difference from the pure spherical wave. A 
suitable correction here is the split Hanning function 



with Xj^ and y^ suitably chosen transition values. This 
treatment is characteristic of many numerical operations 
(such as digital filtering) involving discrete Fourier 
transforms . 

In the calculations reported here, we have generally 
employed a 64 by 64 grid. At the point where the plume 
subtends its maximum angle, between on-third and one-half 
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of the grid dimension lies within the plume. Thus, 20 to 30 
radial points are given within the flow field (less at 
other stations) . This allows a flow field resolution 
comparable to that of a normal flow field calculation. In 
some cases we have experimented with a 128 by 128 point 
grid; the solution for the mean plume is only slightly 
different from that evaluated with the smaller grid when a 
mean (time- aver aged) flow field is used. A 32 by 32 grid, 
however, usually causes some loss in quality. 
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5.0 EXHAUST PLUME FLOW FIELDS 

The flow fields for both the Shuttle ax3,d Titan III-C 
exhaust plumes were furnished by Lockheed Missii.es 
Space Company, Huntsville Research and Engineering Center. 

Two separate codes are used to evaluate conditions in the 
portion of the plume which is expanding to ambient pressure 
and in the downstream pressure-equilibrated trail. The code 
used in the upstream region (RAMP; Penney, et al., 1976) is 
essentially a viscous method-of-characteristics solution; the 
code used downstream (LAMP) is an extended version of the LAPP 
code developed at Aerochem Research Laboratories (Mikatarian, 
et al., 1972). Both codes employ a two-equation turbulence 
model to describe the (time-average) mixing of the exhaust 
with air. Secondary combustion and ionization-deionization 
chemistry are included in both codes. 

The major limitations to these calculations are, 
first, the assimiption of axial symmetry, and second, their 
treatment of a mean rather than an instantaneous flow field. 
The large separation of the SRM's on both Titan III-C and 
Space Shuttle gives their exhaust plumes a pronounced three- 
dimensional character, especially at lower altitudes. Figures 
17 and 18 illustrate the interaction between the exhuast 
plumes at two different altitudes. At low altitude the 
plumes are essentially separate. Such interaction as does 
occur is primarily due to mixing, and in regions in which 
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Figure 18. Estimate of SRB plume boundary locations at 40 km. 



the concentration of exhaust gas is low. Under these 
circumstances it is reasonable to evaluate the two plumes 
independently and to calculate the electrical properties 
in the overlap region by linear superposition of the two 
flow fields. At high altitude {near SRM burnout) the 
plume size is large compared to the motor separation, and 
it is convenient to treat the combined plume as if it issued 
from a single nozzle. This treatment has been shown to 
give good overall plume shapes at very high altitudes, 
although, of course, the details of the region near the 
vehicle are not properly described. At intermediate altitude 
neither treatment is very satisfactory. However, since 
low-cost procedures for treating three-dimensional plumes 
including mixing and chemistry are not available, we have 
adopted the expedient of running the propagation calculation 
for both superposed and merged plume flow fields, in the 
expectation that the difference between them will provide 
some indication of the reliability of the result. 

The plume calculation predicts a mean flow field which 
IS the ensemble (or time) average over a large sample of 
instantaneous flow fields created by the somewhat random 
processes of turbulent convection. The pattern of these 
instantaneous flow fields is irregular on the scale of a 
mixing layer thickness as a result of the large eddy structure 
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of the flow {Fxgure 11) , although when averaged the flow 

properties are smoothly varying on this scale. The trans- 

I 

mitted signal passes through a large eddy in a time ^/c , 

where i is an eddy dimension. A rough estimate for the 

overturning time of this eddy is 2itA/Au , where Au is 

the velocity difference from one side of the eddy to another 

(of order 1 km/sec) . Thus, the signal moves through (or past) 

-5 

the eddy in a fraction of order ttAu/c 10 of its characteristic 
rotation time. The signal thus sees a succession of instan- 
taneous flow structures; plume-induced noise is essentially 
the result of differences in signal properties resulting 
from passage through and around these structures. In 
addition, the mean transmitted signal, being a nonlinear 
function of the flow structure, is not the same when treated 
as the result of interaction with a time-averaged flow 
field as it is when (correctly) treated as the average of 
many interactions with different instantaneous flow fields. 

Because the propagation calculation is relatively 
efficient, it becomes feasible to attempt to model this 
process by constructing several realizations of the turbu- 
lent plume, propagating the signal through each of them, 
and then (incoherently) averaging the results. If we could 
simulate the turbulent field exactly, and perform the propa- 
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gatxon calculation on a fxne enough mesh, this would provide 
an exact solution for the time-dependent signal strength, 
from which we could then infer not only, mean signal level, 
but also many of the noxse characteristics induced by the 
plume. This approach has proven quite fruitful in investx- 
gatxons of the propagation of optical and RF beams through 
the atmosphere. In the present case, we are limited in 
the resolution available , and we can only hope to model 
some of the grosser features of the effects of turbulence. 

Although gross, these features are of interest to 
signal propagation. Inhomogeneities in the plume are ex- 
pected to have a measurable effect on the signal when the 
scale of the xnhomogeneity , 5 ,, times the difference between 
the wave vectors of the scattered and incident signal, Ak, 
xs of order 2,ir. For small-angle scattering Ak kA0, where 
A0 xs the scattering angle, and so Jl k A0/2 it = Z A0/X 1. 

For the range safety signal \ = 0.74 meters, and for S-band 
X — 0.14 meters. Accordxng to field measurements, scattering 
angles from 0 to about 0.25 radians are of interest. The 
scales of xnterest are therefore greater than about 3 meters 

* 

To the order of the approximations invoked in developing 

the propagatxon algorithm. 
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for the Range Safety signal and greater than about 0.5 meters 
for the S-band link. Scales of order 1-2 meters are typi- 
cally resolvable in the propagation calculation near the 
region where the pliime subtends its greatest angle with 
respect to the roll axis. We should therefore expect that 
we can obtain a reasonable estimate of the effects of 
turbulence on the signal at VHP or UHP frequencies, and a 
somewhat cruder estimate at S-band’. Some of the noise 
features should also come from such a calculation, although 
actual simulation of the time-dependent signal strength would 
require a great deal of computation. 

Because the pl^lme is overdense and highly absorbing 
except in the outer regions of the mixing layer, we should 
expect some differences between turbulent effects in plume 
interference and those in (for example) optical propagation 
through the atmosphere. A prominent feature of this region 
is convection of material from the inner regions of the 
plume by large eddies (of the order of the mixing layer width 
in size) . The flow here is essentially intermittent, so 
that time-resolved measurements of local properties at some 
fixed position would show periods of high electron density 
alternating with periods in which the electron density was 
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very low. The probability distribution function of the 
electron density would tend to be bimodal, being the sum of 
an ambient distribution (essentially a delta function) and 
a function characteristic of regions further inside the 
mixing layer. This bimodality of the distribution function 
IS clearly shown in measurements of temperature fluctuations 
in a mixing region between hot and cold air by La Rue and 
Libby (1974) . 

It has often been observed that the boundaries of 

turbulent shear flows exhibit a sharp interface between 

regions of turbulent and non-turbulent flow. The properties 

of this interface have been the object of considerable 

study; La Rue and Libby (1976) have recently summarized 

and extended these studies. Here we can consider the small- 

scale turbulence or some other property of the turbulent 

fluid (elevated temperature in the case of La Rue and Libby) 

as a marker for fluid on the inside of the interface. It 

IS found that many of the statistical properties of the inter- 

* 

face (including the fact that the intermit tency factor is 
an error function of distance transverse to the mean flow, 
that the probability distribution of the crossing frequency 
at a point is essentially Gaussian, and the probability 
distributions of the durations of turbulent and non-turbulent 

* 

Defined as the fraction of time that fluid at a given 

position is turbulent or hot. 
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fluid passage) are well described by considering the position 

of the interface x^. to be a Gaussian-distributed random 

variable. (Alternatively, we could consider the displace- 

ment from its mean position to be the Gaussian variable.) 

Although certain fine details of the interface are not so 

* 

well described by this model, it is an excellent approxi- 
mation for the gross features. 

These observations suggest that a productive treat- 
ment of the turbulent structure of the plume could be 
developed by treating the electron density and collision 
frequency as fluid marker properties which are convected 
about under the influence of a turbulent field in which 
the displacement from a mean position is a Guassian random 
variable. According to one form of the central limit 
theorem, we can produce a Gaussian random variable by con- 
structing a set of realizations in physical space based on 
a Fourier representation whose power spectrum is well- 
behaved and whose phases are a set of random numbers uni- 
formly distributed between 0 and 2 tt. We have adopted this 
approach. 

The turbulent field displaces fluid with certain 
assigned properties from the exoectation value of its position 
<x> to a displaced position <x> + D(x,t) . The turbulent 

ii 

Such as interface slope, which is important to establishing 
the mechanism of entrainment. 
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displacement D is the path integral of the turbulent 
* 

velocity u(x,t) : 


D = 



( 5 - 1 ) 


We can write this in terms of the time derivative of D as 


u 


dD 

dt 


( 5 - 2 ) 


where the derivative is a substantive derivative (following 
fluid elements) : + u*V . 

The remainder of this section is concerned with 
developing a model for the turbulent displacement D . In 
the interests of computational tractability, a number of 
simplifying approximations have been introduced. Most of 
these are justified by approximate physical arguments rather 
than by rigorous derivations . We have attempted to retain 
certain features which could be important, including the 
possibility of correlations between displacements which exist 
at different axial locations in the plume. We do this by 
considering the power spectrum of the velocity and dis- 
placement fields to evolve as fluid is convected downstream 
by the mean flow, and by assigning a set of frequencies to 
the Fourier modes representing the displacement. 
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The evolution of the spectrum presents an immediate 

problem. In principle, we could set up and solve a transport 

* 

equation for the spectrum as a function of position. In 
practice this turns out to be difficult, and we adopt 
instead a procedure in which the local spectrum is deter- 
mined by some overall fluid properties at a given value of 
axial distance Z. On the basis that the turbulent region 
in the plxame is somewhat like a self-preserving mixing 
layer between fluids of different velocities, we expect the 
turbulent velocity spectrum to evolve with distance down- 
stream as follows; Shortly downstream from the origin, 
the mixing region dimension is small, so that the lowest 
wavenumber k^ of the velocity or displacement spectrum is 
large. (This wavenumber is inversely proportional to the 
outer or integral scale of the turbulence; the relation is 
derived in Appendix II.) As the mixing layer grows, the 
cutoff wavenximber decreases; since the maximum scale increases 
linearly with z , we should expect to find k^ decreasing 
as z However, if the mixing region is self-preserving, 

the turbulent intensity must be the same at all stations, 
so that the area under the curve representing the turbulent 
spectrum must be constant. Thus, the spectrum should evolve 
* 

Here we would consider the power spectrum to represent 
spatial and temporal scales small with respect to the 
distance over which the spectrum changes appreciably. 
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as shown in Figure 19. Close to the origin the spectrum 

is made up of small, highly energetic eddies, while farther 

downstream the eddies become progressively larger, but less 

energetic. The turbulent displacement cannot exceed a 

scale proportional to the mixing layer dimension, so that 

there must exist a minimum wavenumber at which the 

o 

spectrum either cuts off or saturates . 

Each spatial Fourier component of the displacement 

will have a correlation time inversely proportional to 

Its frequency, or an axial correlation length proportional 

to f where U is the downstream convective velocity. 

Strong correlations should exist for AZ 1 • while 

essentially no correlations should exist for AZ Wj^/U >> 1 . 

We should expect some effects on the signal if the turbulent 

structure at successive phase screens is correlated. To 

allow for this, we assign a frequency to each of the Fourier 

modes and allow the phases of the modes to change as Wj^t . 

The evolution of the displacement with axial distance 

is evaluated as follows: Suppose that all axial (i.e., z- 

direction) velocities are equal at every point on an x-y 

plane and that they are much larger than all perpendicular 

velocities. Under these conditions we can define a time 

t' measured by an observer moving downstream with the flow, 
3 3 3 

such that. Kirr ~ ^ also assume that we can 

o t o "C 0 Z 
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neglect the convective terms involving the horizontal 
velocities (a reasonable assumption if the mean axial 
velocity is large with respect to lateral turbulent 
velocities) , then we can replace the total derivative in 
Equation (5-2) by a partial derivative, and write 


u- 



(5-3) 


The space-time Fourier transform of this is 








(5-4) 


where is the frequency of the displacement Fourier 

component . 

We now return to the concept that an eddy overturning 
time (its period, x ) is its circumference divided by its 
velocity, or 


T SIS ttX/u 


(5-5) 


In terms of wavenumber and frequency. 


oj «=«= k Uj^ (air) 
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where a is a constant of order unity. Thus the displace- 
ment Fourier components are 

-1 

k (atr) k (5-6) 

where k is a unit vector in k-space. In practice, we 
construct D(x) as a discrete Fourier series. The squared 
amplitudes in these series represent the contribution to 
the total power of the displacement field from a range of 
wavenumbers between k and k + Ak . Thus the discrete 
values of the amplitudes of are 

IDj^! =k (ar) k'^ (Ak^Aky)^ (5-7) 

where Ak and Ak are determined by the grid dimensions . 

We now have an expression for the Fourier components of the 

turbulent displacement and an expression for their frequency. 

Using these expressions we can construct a realization of 

the turbulent displacement field and advance it in time 

(or alternatively in z) as it moves down the plume. 

We expect the three-dimensional turbulent velocity 

-11/3 

spectrum ^^(k) to obey the Kolmogorov k law in 

the inertial range and to have a low-wavenumber cutoff 
or saturation corresponding to the largest eddy size (of 
order of the lateral dimensions of the turbulent region) ; 
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the dissipation scale is well below any scale we must consider. 
A reasonable expression for the turbulent velocity spectrum 
is 


$ (k) = c(k^ + k^) 
u o 


(5-8) 


where the long-wave saturation wavenumber k^ is expected 

to be a steadily decreasing function of time. The constant 

c is chosen by requiring that the mean square tmrbulent 
2 

velocity <u* > be obtained when the spectrum is integrated 
over all wavenumbers. Thus 


c 


<u'^>/ 


I 


f (k) dk 
u 


or explicitly. 


r(ll/6) k^^^ <u'^> 
2tt r(3/2) r{l/3) 


0.512 k^/^ <u'^> 
o 


(5-9) 


(5-10) 


A computational procedure for evaluating realizations 
of the turbulent flow field compatible with the arguments 
given above is the following; At the beginning of each 
separate propagation calculation, a realization of the tur- 
bulent displacement field in Fourier space is calculated 
by assigning amplitudes, phases, and frequencies to a set 
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of modes. We assign N modes to each wavenumber k ; each 
of these N modes has a different assigned frequency. The 
mode amplitudes of the x components of the displacement are 

D°(k^ , ky , »^.) = k^[Ak^ iky]'^ a,r(Nk3)-^ 

X exp(2ir i R) (5-11) 

where R is a random ntamber uniform in the interval (0,1) . 
The y components are evaluated from a similar expression. 
The frequencies of these modes are evaluated from the 
velocity spectrum as follows: We first calculate an expec- 

tation value of the velocity at a given wavenimnber as 

<u,> = (k) k Ak Ak (5-12) 

u X y 

» 

where is the three-dimensional velocity spectrum. 

Equation (5-8) . The expectation value of the frequency is 
then 


<w^> = k<u^>/aiT . (5-13) 

The frequencies of the N individual modes at a given 
wavenumber are assigned by sampling from a Gaussian distri- 
bution. Explicitly,, the relation used is 
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(5-14) 


^ - <Uj,> [1 + R] 

where R is a Gauss ian-distributed random number having 

zero mean and unit variance. 

We have not yet taken account of the long-wavelength 

saturatxon of the displacement spectrum. The saturation 

wavenumber is inversely proportional to the integral 

scale of turbulence, , which we expect to be roughly 

proportional to the width of the mixing region. From 

explicit calculations using a set of turbulence model 

equations (supplied by Lockheed) , we deduce that 

« 0.05 L . , where L . is the width of the mixing 

T mix mix 

layer. Thus at each station we can set 


k 


o 


0.75/L^ 


(5-15) 


This results in a scale for the large eddies (27r/k^) 
of the order of 50% of the mixing layer width, which appears 
reasonable. The saturation factor is simply applied to 
the spectrum at each station; in addition, a phase factor 
containing the frequency is also applied. Thus, 


°x^^x ' ^y ' “ik '■ “^x^^x ' ^y ' “ik^ 


i Z/U 


[k^(z)+k^]^/^ 


(5-16) 
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and simxlarly for D^, . The total displacement component 
at a given value of k is the sum over the N modes: 

i 

The resulting Fourier amplitudes are then subjected 
to an inverse fast Fourier transform which provides the 
displacements as functions of position in physical space. 

The electron density and collision frequency at a point (x,y) 
are assigned the values at the displaced positions 
(x-D , y-D ) in the reference flow field calculation. This 
reference flow field should be chosen such that the expec- 
tation values of the electron density and collision frequency, 
averaged over an ensemble of realizations of the turbulent 
displacement field, are the local values predicted by the 
mean flow calculation. In general, this will mean that the 
reference flow field will be different from the computed 
mean flow: only in special circumstances will the average 

over all turbulent displacements be equal to the local value. 

If the probability that the turbulent dxsplacement is 
between D and D + d5 is P ($) dB , then the mean 
electron density is 

<N^(x)> = (x-B) P(B) d^B (5-18) 
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(x) ->■ 

where (x) is the electron density of the unknown 

reference flow and where we have assumed that P(6) is nor- 
malized to unity. Since this has the form of a two-dimensional 
convolution, we can write the Fourier transform of as 

<N {ic)> = (ic) P(ic) (5-19) 

6 0 

and thus the Fourier transform of the reference field is 

(it) = <N^ (it)>/P(ic) . (5-20) 

Thus if we know P (D) , we can in principle find the 
reference electron density field by transforming the given 
mean field, dividing by the transform of P (x) , and then 
performing an inverse transform on the result. 

It is known that such deconvolutions are often un- 
successful as the result of numerical noise, which results 
in the introduction of spuriously large values at some 
wavenumbers for the reference field (those where P($) 
is small). In addition, test calculations indicate that 

the convolution has only a small effect upon electron density 

2 

profiles unless <D > is a significant fraction of the 
profile width. For these reasons we have neglected the 
convolution error in our turbulence calculation, with the 
result that our reference field — the mean field — is somewhat 
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more extensive than it should be. Matched filter decon- 
volution techniques are capable in principle of compen- 
sating for noise errors, at the cost of some errors in 
localization. These remain a subject for investigation in 
future work. 

Because the grid of the propagation changes in scale 
as the calculation proceeds as a result of applying Talanov's 
tranformation, we carry out the calculation of D(x) in 
physical space on a constant reference grid. This grid is 
chosen to be identical to the computational grid at the 
station where the plume subtends its greatest angle with 
respect to the roll axis. Values of the electron density 
and collision frequency are first calculated on this reference 
grid and then interpolated to the grid of the propagation 

N 

calculation. Our calculation evaluates a displacement 
everywhere in this reference grid, including the potential 
flow outside the plume. We could avoid this either by 
imposing a damping factor on the displacement at some 
distance from the mixing region or (possibly) by intro- 
ducing a phase coherence in D{k) which would limit D(x) 
to the mixing region by phase cancellation elsewhere. On 
the other hand, the electron density outside the plume is 
essentially zero, and convection of it produces no net 
effect unless it is convected into the plume. We therefore 
simply allow D(x) to be evaluated on the entire grid. 
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An example of the effects introduced by this procedure 
is shown in Figure 2 0. The contour plots on the top show 
the imaginary part of the refractive index for the mean flow 
(in this case for two superimposed plumes) at ^a gxven axial 
station; the plots on the bottom show the xndex contours 
resulting from applying a turbulent displacement. The re- 
sultxng irregularity is clearly apparent, and appears to be 
of reasonable size in relation to the overall plume size. 

Note the sharpness of the edge of the plume as evxdenced 
by the close spacing of the lower contours (two plots at 
rxght) . 
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Figure 20. Contours of the imaginary part of the re- 
fractive xndex in turbulent and non-turbulent 
exhuast plumes (turbulent beneath) of a 
Titan III-C at 20 km altitude. Left plots 
show higher contour levels. In each plor, 
dashed contours are lower. 






6.0 DEMONSTRATION AND COMPARISON WITH EXPERIMENT 

We have applied these procedures to the exhaust plumes 
of the Titan III-C and the Space Shuttle. In each case we 
calculate the signal pattern from the antenna to the ground, 
and employ reciprocity to identify the signal strength as 
that resulting at the antenna due to transmission from a 
ground station. Because the near field antenna pattern is 
not available, we have assumed that the antenna is an iso- 
tropic radiator. No blocking by the vehicle is assumed 
except in the plane of the vehicle base. The signals that 
we calculate on the side of the vehicle away from the 
antenna are therefore not to be trusted. 

Because some attenuation data are available for the 
Titan III-C, we can compare our predictions with experiment 
in that case. The quality of agreement between theory 
and experiment (keeping in mind the condition of the data) 
should provide an indication of the error bounds to be placed 
on our calculations for the Space Shuttle. We also have 
used the Titan III-C as a test case for some of the pro- 
cedures we employ to represent the turbulent, three- 
dimensional exhaust plume using the axially symmetric, 
time-averaged flow field. 
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6.1 Sxgnal Interference Effects 

It will be helpful, and illustrative of several 
features of the plume-signal interference phenomenology, 
to follow the progress of the calculation through an 
example: These calculations are for a Titan III-C plume 

at 36 km. The signal frequency is S-band and the two- 
engine plume is treated as a single merged plimie from an 
equivalent single engine. We begin at a z station near 
that at which the plume subtends the maximum angle as 
seen from the antenna location. Figure 21 shows a computer- 
generated contour plot of the field amplitude |^| in a 
plane normal to a line drawn aft from the antenna parallel 
to the vehicle roll axis (and to the plume axis) . The 
plume axis is not identical with the origin of the plot 
because the antenna is one side of the vehicle (the left 
in this case.) The contours are uniformly spaced on a 
linear scale. The lower three levels are dotted, while the 
upper three are solid lines. A maximum occurs between the 
two complete outer circles, while the "islands" outside them 
are a secondary maximum. Small irregularities in the 
contours are the effects of interpolation of the flow 
field data to a rectangular grid, combined with offset of 
the antenna. The contours at the top are values of 
I'i'I/l’I'ol , where is the field at the nozzle exit. 
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Even at this station, about 10 meters downstream from 
the engine, the field is very strongly attenuated near 
the plume centerline. 

Figure 22 shows a plot of the imaginary part of the 
plume's refractive index in the same plane. Comparison 
of this figure with the previous figure shows that the 
region of the field in which the amplitude is appreciable 
lies mostly outside the plume. The attenuation at the 
lowest contour level is approximately 4 db/meter. The 
result of subjecting the incident field to passage through 
this phase screen is shown in Figure 23, (The screen 
thickness is about 2 meters.) The values of are 

shown for the same levels as in Figure 21, Little has 
happened to the external field, but that at the edges of 
the plume has been absorbed, as evidenced by the greater 
displacement of the low contours from the plume axis. 

The behavior of the calculation at the other stations is 
roughly similar. Figure 24 shows the values of the field 
amplitude just after the last phase screen included in this 
calculation. Some diffraction has evidently occurred, 
since the higher contours are more spread out than before. 
The field strength is lovzer as a result of geometric 
spreading; note the larger values of x and y . The region 
of appreciable phase change (Figure 25) occupies only a 
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Figure 22. Contours of the imaginary part of the refractive index 
(high at left, low at right) at the axial station of 
Figure 21. 
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Figure 23, Contours of constant field amplitude 
after emerging from the phase screen 
represented by the previous figure. 
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Figure 25. Contours of the imaginary part of the 
phase change (attenuation exponent) 
in the last phase screen considered. 





very small portion of the grid at this station; this is a 
consequence of the plume's spreading less rapidly than 
the mesh, which grows in accordance with geometric optics 
under the rules of the Talanov transformation. 

Finally, Figure 26 shows intensity contours 
in a far-field plane (in this case, 50 km 
behind the vehicle) . The intensity is normalized with 
respect to the uniform free-space intensity, and the 
contour levels are in db (volts) . The abscissa and ordinate 
are angles in the x and y direction, in degrees. The angular 
resolution of this calculation is about 1.3 degrees. The 
intensity is given at each of 4096 points in this grid. 

Figures 27 and 28 show the intensity ratio in db along the 
= 0 and ~ ^ axes of the previous figure. Only 

the 6 =0 cut is symmetric; this is a result of the 

initial offset of the plume axis from the nominal propa- 
gation axis in the x direction. The maximum attenuation 
level is about 90 db. 

We have performed a few elementary variations about 
this basic type of calculation illustrating the importance 
of the various elements of the calculation. When the 
electron density ^t every point in the plume is doubled, 
the results are very little different from the nominal 
case. This indicates the relative unimportance of conditions 
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Figure 26. Contours of constant normalized field 
intensity (in db) in the far field 
(Titan III-C at 36 kiti) . 
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Figure 27. Far-field intensity as a function of 
aspect angle 6^ (normal to motor 
plane) for 0=0. 
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Figure 28. Far-field intensity as a function of 


aspect angle 
for 9 = 0 . 
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( in motor plane) 
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within the plume when the absorption is high. It also 
depends upon the relative steepness of ttie electron density- 
gradients near the edge of the plume; if these are sharp, 
then even a siJastantial change in the value of n^ will 
not produce much change in -the position of the plume edge, 
the extent of the shadowed region, and the diffraction 
pattern around it. However, if we increase the dimensions 
of the plume by a modest fraction (say 10-20%) , the 
attenuation pattern in the far field becomes appreciably 
wider and deeper. These trial calculations indicate that 
simple geometric blockage by the plume is the predominant 
effect. 

The irrelevance of internal plume properties to the 
signal is again illustrated by Figure 29. This shows the 
signal attenuation per unit length expected in a Titan III-C 
plume at 36 km altitude as a function of radial distance 
at the axial station shown in Figures 21-23 at P-band and 
S-band frequencies. The dots are the actual points from 
the flow field calculation. Note that although the S-band 
signal is less attenuated than the P-band signal in the 
outer region of the plume, the difference is essentially 
irrelevant. The attenuation levels of interest are less 
than 1 db/meter, and -these low levels are attained only 
at the plume edge. Note also that the effective S-band 
plume at those attenuation levels is not much less wide -than 
the P-band plme. 
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SIGNAL ATTENUATION, db I /meter 




Fxgure 30 shows the geometry of the plume and the 
far-field signal in schematic fashion. The plume is here 
represented by the actual set of phase screens used in the 


calculation (excluding some far downstream which have no 
effect) . Each perpendicular line from the abscissa repre- 
sents a phase screen; the edge is at the 1 db /meter attenuation 
level. The dashed line angled upwards from the antenna 
represents the edge of the shadowed region (umbra) . The 
amplitude of the far-field signal is shown at right, as 
a function of aspect angle. The geometric teirminator passes 
along the edge of the plume for an appreciable distance. 
Multiple scattering at several axial stations undoubtedly 
has some influence on the result. The dominant effects 
of geometry upon the signal are clearly shown. 

Introducing simulated turbulence into this calcula- 
tion has a noticeable effect upon the behavior of the 
plume and signal at intermediate stages, and a less signifi- 
cant effect in the far field. The plume is essentially 
sharp-edged on the grid scale, and the predominant effect 
of turbulence is to distort the edges of the plume. Figure 
31 shows the resulting signal at a station in the plume for 
a specific realization of the turbulent displacement. 

(Compare this figure with Figures 21-23.) The incoming 
field is shown at left, and the outgoing at right. It 
appears that the turbulence at this station has to some 
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Figure 30. Plume-signal interaction geometry. Plot at right shows 
normalized field amplitude as a function of aspect 
angle (S-band) . 
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Figure 31. 


Incident (left) and emergent (right) fields at an axial 
station (cf. Figures 21-23) in a turbulent _ Titan III-C 
plume ( individual realization) at 36 km altitude. 






extent counteracted the effects which occur at earlier 
statxons, and made the field a bit smoother. One might 
expect this behavior if the structure decorrelates quickly 
with axial distance (or convected time) . 

Figure 32 shows the effects of simulated plume 
turbulence on the far-field signal. Here we show several 
plots of field intensity versus aspect angle in the plume 
of the two solid motors for the Titan III-C. (These plots 
would be symmetric in the absence of any turbulent effects; 
the calculation is otherwise identical with that of Figure 
28.) The abscissa is identical with that of Figure 28; the 
ordinate is computer-generated and varies slightly. Peak 
attenuation (essentially tail aspect) varies from about 
80 to 90 db in these calculations. The introduction 
of simulated turbulence has not greatly affected the overall 
signal attenuation pattern, although several details (es- 
pecially the position and shape of fringes) are affected. 

The variability in the predicted signal is greatest at 
tail aspect (standard deviations of 4-6 db) and generally 
decreases toward the edge of the shadowed region. No large 
signal dropouts are predicted by this model of propagation 
in (or better, around) a turbulent medium. 

At lower altitudes the effects of turbulence are 
somewhat greater. Figure 33 shows a realization of the 
incident and emergent wave fields at an axial location 
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Figure 32. Plots of field intensity in db vs. aspect angle Qy 

(in solid motor plane) for several individual realizations 
of simulated plume turbulence. Titan III-C at 36 km. 
Compare with Figure 28. 
















roughly equivalent to that of Figure 31 for the Titan III-C 
at 20.6 km, and Figure 34 shows the far-field amplitude 
variation for a set of realizations. (In these calculations 
we treat the flow field as a superposition of two indepen- 
dent plumes.) The attenuation for the mean flow is shovm 
in Figure 35. Figure 33 indicates that the turbulent 
displacements are somewhat greater with respect to the 
plume outline at low altitude. This is probably a conse- 
quence of the mixing layer's representing a smaller fraction 
of the plume cross section at higher altitude, where a 
greater degree of inviscid expansion takes place. The tail 
aspect signal attenuation in these calculations varies 
between 88 and 102 db. Note that the effects of turbulence 
on the position and intensity of the secondary diffraction 
minima appear to also be greater at lower altitude. These 
effects will have some influence on the AM noise level in 
the signal, although we do not consider this calculation 
sufficiently well checked out to quote a value. The effects 
of simulated turbulence are considerably less at P-band 
than at S-band. 

We have performed a few calculations in which the 
effects of the parameters entering the turbulence simulation 
have been varied in order to determine the sensitivity of 
the propagation calculation to them. The scale of turbu- 
lence (which determines the saturation wavenumber k^ — see 
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Figure 34. 


Plots of field intensity j ip | in db vs. aspect angle _0 
(in solid motor plane) for several individual realizatic 
of simulated plume turbulence. Titan III-C at 20 km. 
Compare with Figure 35. 
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Figure 35. Plot of normalxzed intensity in 
db as a function of aspect angle 
0y (in motor plane) in the far field 
Txtan III-C at 20.6 km. 



Section 5 and Appendix II) is the parameter which has the 
greatest effect. Doubling or quadrupling the scale appears 
to introduce additional noise and 4-6 db greater signal 
attenuation at all angles where it exits. This finding 
should be taken as a caveat on our results. The scale in 
the outer regions of the mixing layer could be somewhat 
different from what we have assumed. It may he possible 
to determine effective values of the parameters by comparison 
with various data. Because coding of the turbulence model 
was not completed until late in the study, we have not 
done so. 

6 . 2 Titan III-C Calculation and Comparison with Experiment 

Calculations for the Titan III-C plume were performed 

at 20.6 and 36.0 km altitude, corresponding to times 80 and 

110 seconds after liftoff, in order to compare the model 

predictions with field data. The first time corresponds 

approximately to the onset time for appreciable attenuation 

in Poehler ' s measurements. Only the solid motors (Stage 0) 

are burning at this time. The Stage I engines burn a 

mixture of hydrazine and unsymmetrical dimethylhydrazine 

(UDMH) with nitrogen tetroxide. Unlike the hydrogen- 

oxygen fuel used in the Shuttle orbiter stage, the Titan III-C 

* 

Stage I propellant produces a somewhat ionized exhaust 

and its plume produces attenuation at higher altitudes 

where the SRM's have been dropped (see Figure 3) , and we 
_ 

H.A. Poehler ascribes the observed attenuation to the release 
of alkali metal impurities from the ablative shirt of the 
nozzles . 
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wish to avoid this complication. We thus do not treat the 
pliime after Stage I ignition. 

Figures 36-39 show calculated attenuations for the 
smooth flow field of the Titan III-C plumes at these 
altitudes. At the lower altitudes we have treated the 
plume as the superposition of two independent plumes from 
the two engines. At the higher altitude we have treated the 
plume as a single "merged" plume, i.e., as if it came from 
a single nozzle with the combined thrust of the two motors, 
by simply multiplying all plume dimensions by a factor of 
/2 . In doing so we make errors in plume chemistry and 
local electron densities inside the plume; our previous 
calculations indicate that this is permissible. 

These calculations are compared with Poehler's measure- 
ments in Figures 40-43. Since the observations are made at 
a variety of azimuthal angles, we show calculations both versus 
0^ (normal to motor plane) and 6^ (in motor plane) . The 
agreement between the calculation and the experiment is generally 
satisfactory at P-band and less so at S-band. One feature 
of the calculations which does not appear in the data is 
the strong attenuation of the S-band signal near tail aspect. 
Since the calculation is in fair agreement with measurements 
at angles greater than 5°, we believe that this strong tail- 
aspect attenuation is real. The data do not show it because 
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Figure 37. S-band attenuation vs 
Titan III-C at 36 km. 
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Figure 41. Comparisons between predictions and 
field data: Titan III-C at 36 km, 

S-band. 
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Comparisons between predictions and 
field measurements: Titan III-C at 
36 km, P-band. 


the fixed location of the ground stations and the launch 
geometry did not allow any measurements to be made at angles 
less than about 5 degrees. 

Our investigations of the effects of simulated 
turbulence in the previous section indicated that large 
effects were not to be expected for mean signal attenuation. 

An average over a modest ensemble (10 realizations) of 
the turbulent displacement field did not produce any great 
difference from the smooth- flow calculations, although 
some smearing of the interference fringe effects near 
10-12” at S-band was observed. 

These comparisons indicate that our model has well 
reproduced the phenomenology of RF interference by large 
solid-motor plumes. The agreement with data is generally 
acceptable, although the high-altitude S-band predictions 
apppear to fall substantially below the measurements. Because 
our numerical resolution is not sufficient to resolve turbu- 
lent structures of the sort which would scatter S-band radiation 
strongly through the angles involved in bending the plume- 
skimming ray back toward tail aspect, it is conceivable that 
subgrid-level scattering processes may be involved in producing 
the high measured signals under these conditions; this subject 
remains to be explored. In any case, the calculation appears 
conservative, and it can be applied to the prediction of 
interference effects of the Space Shuttle plumes with reasonable 
confidence . 
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7.0 RESULTS FOR SPACE SHUTTLE 


7 .1 Results 

We have performed several calculations of Shuttle 
SRB plume attenuation at the Range Safety signal frequency 
of 416.5 MHz. Our results are for the antennas on the 
SRB's; sxmilar effects are expected for the external tank 
(ET) location. These calculations specify field intensxties 
on a 64x64 grid in a far-field plane normal to the vehicle 
roll axxs. Detailed results will be furnished separately 
to NASA/MSFC; xn this report we show only selected features. 

Figures 44-46 show attenuation results as functxons 
of aspect angles 0^ and 9^ . Here the 6 ^^ aspect is 
in the plane normal to the SRB ' s and the 0^ aspect is 
in the plane of the SRB's. The axis is a line parallel 
to the roll axis passing through the antenna- Because the 
antenna is offset in both x and y from the roll axis, 
no plane of symmetry exxsts in these calculations. The plume 
xs treated as a linear superposxtion of two independent 
SRB plumes at the two lower altitudes , and as a single 
"merged" plume at the higher altitude. These calculations 
are for the mean flow field. Because of vehicle blockage 
effects not included in our calculations, only results on 
the sxde on which the antenna is located are likely to be 
meaningful. This is the side of negative 9 ^ and 0 y . 
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Figure 46. Predicted Space Shut 
as a function of the 
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Obtaining reliable results on the far side of the vehicle 
would requite dropping our assumption that the antenna 
pattern is isotropic and (most probably) embarking on 
detailed analysis of vehicle skin currents / a task outside 
the scope of this study. 

In addition to these calculations for the nominal 
flow fields, we have examined the effects of turbulence and 
of variations in descriptions of turbulent mixing and 
electron chemistry in the edge region. In the case of the 
latter two quantities, Lockheed has estimated a correction 
to the effective radial position of the plume boundary 
which depends on axial distance and (in the case of mixing) 
on turbulence. The effect in both cases is to displace 
the boundary slightly inward, i.e., toward the plume 
centerline, from the nominally calculated value. The change 
in mixing model makes electron-rich material diffuse less 
rapidly, while the changes in chemistry result in more 
rapid electron loss in the outer regions of the plume be- 
cause of attachment to atmospheric species. 

Figures 47-49 show contour plots of the normalized 

2 2 2 
far-field Intensity, j j /|^q 1 , m db where 

IS the intensity in the absence of a plume, for the 

nominal plume calculation. The lower left quadrant of 

these and succeeding plots corresponds to angles in the 

quadrant of the antenna with respect to the roll axis. 
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Figure 47. Contours of normalized far-field 

xntensity for Space Shuttle at 9.3 km 
altitude. Nominal mean flow field - 
Contour levels on this and succeeding 
plots are in db . 
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Figure 48, Contours of normalized far-field 

intensity for Space Shuttle at 25 km 
altitude. Nominal mean flow field. 
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Figures 50-52 show similar contours for calculations in 
which Lockheed's recommended plume boundary position 
variation has been applied to the nominal flow fields. 

In these cases the plume with variations in mixing rate 
is on the left, and the plume with variations in edge 
electron chemistry on the right. (The easiest way to 
compare these with the nominal contours is to make copies 
and overlay the contour plots, since the contour levels are 
generally similar.) In general, applying these variations 
to the flow field results in a slight decrease in attenuation 
(or increase in the transmitted signal) . Local attenuation 
increases can occur as the result of changes in the 
position of interference fringes. Since the effect of tur- 
bulence is to smear out these fringe patterns, we do not 
believe that these modest local increases in attenuation 
should be considered important. 

The effects of turbulence on the normalized far- 
field intensity are shown in Figures 53-55. These calcula- 
tions are the average amplitude over ten realizations of the 
turbulent plume. Including simulated turbulence causes a 
slight broadening and deepening of the attenuation pattern, 
similar to the effects shown in Section 6.1 The effects 
are greatest at tail aspect and somewhat greater at the 
higher altitudes. In general, large effects are found only 
when the signal loss due to the mean flow is high (S70 db) . 
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Figure 50. 


Contours of normalized far-field intensity for Space Shuttle 
at 25 km altitude. Left, mixing variation; right, edge 
kinetics variation. 
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Figure 51. Contours of normalized far-field intensity for Space 
Shuttle at 25 km altitude. Left, mixing variation; 
right, edge kinetics variation. 
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Figure 52, Contours of normalized far-field intensity for 

Space Shuttle at 37.5 km altitude. Left, mixing 
variation; right, edge kinetics variation. 
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Figure 53. Contours of normalxzed far-field 
intensxty for Space Shuttle at 
9.3 km altitude. Average of 10 
realxzations of simulated turbulence. 
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Figure 54. Contours of normalized far-field 

intensity for Space Shuttle at 25 km 
altitude. Average of 10 realizations 
of simulated turbulence. 
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1.2 Comments and Conclusions 


Comparisons of our signal interference calculations 

for Space Shuttle during SRB burn with those for Titan III-C 
th 

0^ stage indicate that, at any given altitude, the Space 

Shuttle plume induces a greater signal loss than the 

Titan III-C plume at UHF-VHF frequencies. This is a result 

of the generally increased geometric blockage for the 

Shuttle, as evidenced in Figures 8 and 9. Although we have 

not performed calculations at S-band, we believe that the 

interference effects for Space Shuttle will also be greater 

at high frequencies. Signal losses of 30 db or greater 

have been said to cause difficulty with the Range Safety 
* 

system. The aspect angles at which these attenuation 
levels are predicted are about 1°, 9°, and 15-20°, depending 
on azimuthal angle, for altitudes of 9.3 km, 25 km, and 
37.5 km, respectively. The implications on location of 
transmission sites depend upon whether the system must 

transmit a signal under nominal or "worst-case" trajectory 
conditions. The range of excluded aspect angles 

increases with increasing altitude and should continue 
to do so until SRB burnout. A higher or lower value of the 
maximum acceptable signal loss would result in onset of the 


Lee Malone, NASA/MSFC, private communication. 
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attenuation problem at a lower or higher altitude, respec- 
tively, and a greater or lesser range of excluded aspect 
angles at the higher altitudes. 

In applying these results to evaluating the locations 
of alternate transmission sites for Range Safety transmitters 
(to insure that at least one transmitter has a low-loss 
path to the vehicle) , a reasonable degree of conservatism 
would appear to result in preferring the predictions based 
on the nominal flow field over those in which the mixing rate 
and edge Kinetics variations have been included. The com- 
parisons of the Titan III-C predictions with measurements 
indicate that the mean flow predictions are close to 
measured values. The effects of simulated turbulence cause 
a modest broadening and deepening of the attenuation 
pattern. Because coding of the turbulence model was com- 
pleted late in the study period, we have not had much time 
to study its effects on the signal under various parameter 
variations. It is also desirable to compare the results 
with a greater variety of data, such as those to be obtained 
in static test experiments. We therefore recommend that 
our mean-flow results be used in design analyses, but that 
a 5-db safety factor be applied to any calculated result in 
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which the nominal predicted signal loss is less than 30 db. 
This procedure not only is consistent with the effects of 
including simulated turbulence, but also would result in 
our calculated results for Titan III-C lying at the lower 
edge of the scatter in Poehlers ' data. 

In all our calculations we have assumed that the 
antenna pattern is isotropic. We have already noted that 
our results are probably unreliable on the occluded side 
of the vehicle. If the near— field pattern were known in 
the vicinity of the vehicle, we could have used it as an 
initial condition to the propagation calculation; however, 
only the far-field (Fraunhofer) pattern is known. 

The total signal loss is a combination of antenna 
pattern losses and plume-induced losses. Because the trans- 
mitted signal is diffracted around the edges of the plume 
rather than propagated directly through it, the resultant 
signal loss at some point in the plume shadow is not simply 
the db sum of the antenna loss and the plume loss at that 
point. Calculating the signal loss by such a procedure 
would be overly conservative in the case that a region of 
high antenna loss coincides with a region of high plume 
loss. One could conceive also of conditions under which it 
IS optimistic. A more realistic procedure is to add the db 
plume loss to the db antenna loss at the plume edge. Because 
the signal near small aspect angles, where high plume-induced 


123 



losses occur, xs the contribution of rays diffracted over 
a region of the plume edge, it is reasonable to calculate 
the effective antenna contribution as a linear (not 
logarithmic) average over some range of azimuthal angles 
at aspect angles where the plume- induced loss is 3 db. 

The range of azimuthal angles will depend upon the shape 
of the plume pattern and the location of the point in 
question. For small aspect angles, an average over 90° 
in the quadrant about the roll axis containing the antenna 
appears reasonable. 
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APPENDIX I: The Parabolxc Wave Equation and the 

Split Operator Solution 

It IS easy to show, starting with Maxwell’s equations, 
that the basic wave equation for the electric field ampli- 
tude is 


1 

2 

c 




(l-l) 


where e is the dielectric constant. Letting the time- 
dependent part of the quasi-steady solution be given as 
e^^^ where w = ck , we obtain Helmholtz ' equation 


2 

E + 


k e E = 0 


( 1 - 2 ) 


Nov 7 assume that propagation is predominantly in the direction 

of the z axis and that the field is represented as the 

X Ic z 

product of a sinusoidal terra e and a slowly varying 
terra >j7(x,y,z); then 

ijj] + [V^ + k^ £ ] e^^^ = 0 . (1-3) 

Carrying out the z differentiation and neglecting the 
second derivative of the slowly varying quantity rp gives 

2i k 1^ = ip + k^(e-l) = 0 (1-4) 

o Z X 
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which is the parabolic wave equation. 

A formal solution to this equation over a propa- 

^ ^ n n+1 

gation step z z is 




n+1 


exp i 


1 

2k 


2 2 
Az + k 


/ 


n+1 


(e-1) dz 


(1-5) 


n 


For any analytic function (e-1) such that a Taylor 
series in z exists, this is equivalent to second order in 
z to the symmetrized split operator form 


-^AzV^ ik. 7 :pr- i . „2 

,n+l 4k X - ^Az(e-l) - jjjAzV 

V = e e e ip 


( 1 - 6 ) 


The proof below is taken from Fleck, Morris, and Feit (1975) 
Consider the somewhat simpler equation 


II = a(z) = (V^ + b) . 


(1-7) 


Formally, 


r ' z 


Ip (z) = exp 


I 


a (z ' ) dz ' 

to 


( 1 - 8 ) 


where the exponent is taken to mean 


exp 


^ z 

I 


a(z') dz' 




z z z ' 

= 1 + J" a(z')dz' + J" a(z')dz' f a(z")dz" 
o o o 

+ (1-9) 
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Substituting for a{z) and keeping terms through second 
order in z we have 


r 2 


exp J a ( z ' ) dz ' 


= 1 + V 


* J 


b (z) dz ' + V 


2 zf 
1 2 


z ' 


+ V' 


Jdz. I 


£d 

I 


b(z")dz" + / b{z')z’ vj dz' 


z z ' 

+ J" b(z')dz' y* b(z")dz" 
o o 


Also to second order in z , 

rZ, 


1,2^ fb(zMdz' 1^2 


2 1 o 
e e 


2 


= 1 + V 


/ 


b(z’ )dz' 


z, z. z 

J b(z')dz' J b(z")dz" + I J* b{z')dz' 


I / 


) dz' 


and therefore 


(I-IO) 


(I-ll) 
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r z 


exp K 


U 


a(z ' )dz' i ^ 


^V^z I b(z')dz' 


z r ^ 1 

- vj Jdz' |b(z') - j b(z")dz'J 


4i 

- J dz' 


b(z') f - z' vj + ^(z^) . 


(I 


Now if b(z) has a Taylor series representation in z , 

b(z) = b{o) + b' (o) z + b" (o) z^ + (I 


all the remainder terms are third order. Thus, through 
terms of second order in z , the split operator form is 
a solution to the model equation, which differs from the 
parabolic wave equation only by nomenclature. 
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APPENDIX II : Relation of to the Integral Scale Turbulence 

The integral or outer scale of turbulence is defined 
as {see, e.g. , Bradshaw, 1976) 


CO 

^T,x "I / dx/<u^(o)> {II-D 

— Oo 

where <u (o) u (x)^ is the two-point velocity correlation 
function for the x direction. (In isotropic turbulence 
this scale is independent of direction.) Here we have used 
the fact that <u(o) u(x-)^ is a symmetric function of x 
to extend the range of the integral to -<» . A one- 
dimensional power spectrum of the turbulent fluctuations is 
defined as 


CO 

■t'xx'V ■ -t “x'==>> ® 


ik x 

X 


dx 


(II-2) 


so that 




CO 

u^(x) = J 


ik X 

(j) (k ) e ^ dk 

XX X X 


(11-3) 


The three-dimensional power spectrum s (^) is defined as 

2^^ 


»xx<*> = 77^ “x‘*»> ^ 

(27T ) 


-ik»x ,3 -i- 

d X 


(II-4) 
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Here the xx subscript reminds us of the direction of the 
velocities. The one-dimensional spectrum defined above is 
the projection of the three-dimensional spectrum on the 
axis : 


00 CO 
—CO —Co 


4 (k , k , k ) dk dk 
' XX X y z X y 


(II-5) 


(It should be distinguished from another kind -of one- 
dimensional spectrum (j){k) , which results from integrating 
$ (k) over all solid angles at a given value of k = jkl • 
It is this other spectrum ^ (k) which obeys the Kolmogorov 
law in the inertial regime.) In terms of the 1-D and 3-D 
spectra 4> and f we can express L as 


CO oo 

^ - I / 


ik X 


X 


x//l 


^ dk dx//$ „(k) d^ k. (II-6) 


XX 


_oo —oa 


Using the relation (II-5) between the 1-D and 3-D 
spectra, we can write 


CO 00 CO OO 


2 J J J f 


^ ik X 

$ (k) dk dk e ^ dk dx 

XX y z X 


/ 


/' 


—CO —CO —CO —CO 


V ^ 3 ^ 

$ (k) d k 

X 


(II-7) 


Upon changing the order of the integrations, we have 
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L 


T 


CO CO CO 


= 1 T 


///* 

—CO —CO —CO 



CO 


I 


ik X 

e ^ dx dk^ 


—CO 





(II-8) 


The inner integral is 2-n- 6 ' where 6 (x) is the Dirac 

delta function. Thus we can perform the integration 

immediately, obtaining 

oo CO . 

= , r j4.(o , ky , k^) dky dky j*(S) 6 . (II-9) 

— —CO 


If we now assume isotropic turbulence, the spectrum 
IS a, function of k = j ^ I only • Under these conditions 
we can transform to polar coordinates (cylindrical in the 
numerator, spherical in the denominator) and carry out the 
trivial angle integration to arrive at the result 


L = 2tt‘ 


00 

/ 


f (o,K) K d K 


OO 

A/ 


9 (k) k dk 


(II-IO) 


where 


K 


(k^ + 
Y 




(k' 


- k^) ^ 

X 


(II-ll) 


inserting our model expression [Equation (5-8) of the text] 


$(k) 


c 





( 11 - 12 ) 
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we then have 



Thus the saturation wavenumber is approximately the 

inverse of the integral scale of turbulence, an intuitively 
obvious result. 
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